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Abstract 

The Holstein-Hubbard model is examined in the limit of infinite dimensions. 
Conventional folklore states that charge-density-wave (CDW) order is more 
strongly affected by Coulomb repulsion than superconducting order because 
of the pseudopotential effect. We find that both incommensurate CDW and 
superconducting phases are stabilized by the Coulomb repulsion, but, surpris- 
ingly, the commensurate CDW transition temperature is more robust than the 
superconducting transition temperature. This puzzling feature is resolved by 
a detailed analysis of perturbation theory. 
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Introduction. Interacting electronic systems with both electron-electron repulsion and 
phonon-mediated attractive interactions display a rich variety of ground states due to the 
competition between these interactions. In conventional low-temperature superconductors 
(SC), the electron-phonon interaction dominates over the Coulomb repulsion, since the latter 
is reduced from its bare values by the so-called pseudopotential effect (as described by 
Bogoliubov, et. al. Q and Morel and Anderson f2|). The electron-phonon interaction is also 
believed to be responsible for charge- density- wave (CDW) order in weakly doped BaBiC^. 
However, the Coulomb repulsion dominates over the electron-phonon interaction in the high- 
temperature superconducting cuprates, with commensurate spin-density-wave (SDW) order 
occurring in the (undoped) parent compounds. 

Although the Coulomb interaction (in the form of the Hubbard model 03) and the 
electron-phonon interaction (in the form of the Holstein model ||) have both been exten- 
sively studied, there has been only limited work on the combined Holstein-Hubbard model 
H and only one exact theorem for the case of an attractive Coulomb interaction 0. 

The conventional lore for the effect of the Coulomb interaction on a superconductor with 
strong electron-phonon interactions, is that the Coulomb repulsion is reduced from its bare 
values. This is because the electron-phonon interaction is retarded, allowing the electrons 
to attract each other, through the exchange of a virtual phonon, without being at the same 
lattice site at the same time. A quantitative estimate for this so-called pseudopotential effect 
can be made, and the (dimensionless) Coulomb repulsion p(fi)U c is reduced to 

l + p(/i)C/ c ln^ 

where is the electronic density of states (DOS) for an electron (of one spin) at the 
Fermi energy, U c is the bare Coulomb repulsion, W is the electronic bandwidth, and uj£> is 
the Debye frequency. There is, however, no pseudopotential effect for a CDW distortion, 
because retardation effects play a limited role in a static CDW, where the electrons remained 
paired at every other lattice site. Therefore, we expect that the Coulomb repulsion will 
reduce the transition temperatures for CDW order much more than for SC, and that the 



SC phase is thereby stabilized relative to the CDW phase. 

Formalism. The dynamical mean-field theory (MFT) of Metzner and Vollhardt || has 
been employed to exactly solve the Hubbard || and Holstein [HJ models in infinite spatial 



dimensions using the quantum Monte Carlo (QMC) algorithm of Hirsch and Fye ||11|| . The 
Holstein-Hubbard model 

t* 

+ \H^ + M^xl) , (2) 

% 

is written in standard notation: c\ a is a creation operator for an electron localized at site % 
with spin a; := c\ a Ci a is the corresponding number operator; x« is the phonon coordinate 
at site i; and pi is its momentum. The hopping integral is t =: t* /2\fd, the deformation 
potential is g, U c is the Coulomb repulsion, Q is the phonon frequency, and M = 1 is the 
phonon mass. The scaled hopping integral t* determines the energy unit and is set equal to 
one (t* = 1). The effective electron-electron attraction (due to phonon exchange) satisfies 
U = —g 2 /MQ 2 , and competes with the Coulomb repulsion U c > 0. 

The Holstein-Hubbard model is solved by QMC simulation. The QMC algorithm deter- 
mines the local electronic Green's function G(iu n ) at each Fermionic Matsubara frequency 
uj n := 7rT(2n+ 1), by mapping the infinite-dimensional lattice to an impurity problem. This 
mapping has been described in detail elsewhere ||3|,|T2l . The Green's function satisfies 



G(iuj n ) = Fooliun + n - E(iu n )] , (3) 

with S(za; n ) the electronic self energy and F^z) := / dyp(y)/(z — y) the rescaled comple- 
mentary error function of a complex argument [p(y) '■— exp(— y 2 )j\pK is the noninteracting 
DOS]. 

The momentum-dependent susceptibility x(q) = T^Xmn(q) =: rExfli^nn^n) (f° r 
CDW, SDW, or SC order) satisfies a Dyson equation 

Xmn(q) = Xmi^Smn X^il^ mpXpni^) (4) 



for each ordering vector q, with x° the relevant bare susceptibility, and T mp the local irre- 
ducible vertex function. The bare susceptibility for commensurate [q = Q := (71", 7r, 7r, ...)] 
CDW or SDW order is 

x°»(Q):=- . + {iUn L , , (5) 



with a more complicated form for incommensurate wavevectors |13| . On the other hand, the 
uniform bare susceptibility for SC order satisfies 



*g(0) = - ""f , • (6) 

The irreducible vertex functions are extracted directly from the QMC data ||. Figure 1 
displays the lowest-order diagrammatic contributions to the vertex functions in the (a) CDW, 
(b) SC, (c) and SDW channels. 

Transition temperatures are found by starting in the disordered (high-temperature) 
phase, and reducing the temperature until the susceptibility for each ordered phase diverges. 
The highest transition temperature T c determines the initial symmetry of the ordered phase. 

Results. The phase diagrams for the Holstein-Hubbard model with g = 0.5, Q = 0.5, 
U = —1.0, and U c = 0.0, 0.25, 0.5, 0.75, are displayed in Figure 2. These phase diagrams 
are determined by QMC calculations and by a second-order iterated perturbation theory 
(IPT) (there are no detectable phase transitions with the IPT for U c = 0.75). The 
solid dots (lines) depict the commensurate CDW, the open dots (dotted lines) depict the 
incommensurate order, and the open triangles (dot-dashed lines) depict the SC phase for 
the QMC (IPT) calculations (the dashed lines connecting the QMC points are a guide to 
the eye). 

The QMC data displays two types of notable behavior. First, there are no stable in- 
commensurate phases when U c — 0; the incommensurate phases become stable near the 
CDW-SC phase boundary only for U c > 0. The explanation for this is simple: if the SC 
phase was ignored, the CDW phase would suffer a commensurate-incommensurate phase 



transition as T c — > in a similar fashion to the repulsive Hubbard model |15| ; however 



the SC transition temperature is greater than the highest incommensurate CDW transition 
temperature, precluding its appearance at U c = 0. As U c increases, the SC T c drops below 
the maximal incommensurate CDW T c , which allows the incommensurate order to occur. 
Second, even though the region where the CDW phase is stable shrinks as U c is increased, 
the CDW transition temperature is reduced by a smaller factor than the SC transition tem- 
perature, in opposition to the conventional folklore. 

The IPT approximation is reasonable for both the CDW and SC T c 's and for the phase 
boundary between CDW and SC order. The approximate CDW transition temperature is 
again reduced by a smaller factor than the SC transition temperature as U c increases. Thus 
one can understand this effect by studying the weak-coupling formalism. The IPT errs only 
by predicting a large incommensurate CDW-ordered region at U c = which shrinks as U c 
increases, exactly opposite to what the QMC found, and it is unable to predict any finite 
T c 's for U c = 0.75. 

The modification of the CDW transition temperature at half filling is plotted versus U c 
in Figure 3. The CDW phase is followed for < U c < \U\, since a SDW phase is expected 
to be the stable phase for U c > \U\ at half filling (this can be seen at weak coupling 
by comparing the CDW vertex to the SDW vertex in a power series as shown in Figure 1, 
indicating Tsdw > ^cdw when U c > \U\). T C {U C ) is smaller than T C (U C = 0) in the weak- 
coupling regime (g < 0.625) and the curves are nearly linear in U C /\U\, with a decreasing 
slope as g increases. In the strong-coupling regime (g = 1.0), the Coulomb repulsion initially 
enhances the transition temperature (since it reduces the energy of the virtual state formed 
by breaking the bipolaron) before causing a reduction as U c — > \U\. 

Theory. Much of the unexpected and notable behavior found in the QMC and ITP 
results can be illustrated within an analytic approximation. A weak-coupling analysis of the 
CDW and SC transition temperatures is performed in the square-well approximation [where 
the soft cutoffs Q 2 / {Q 2 + {u m — uj n ) 2 } are replaced by hard cutoffs 9(uj c — |u; m |)^(a; c — \ui n \)] 
|TJ|. A first-order calculation is accurate only to the lowest order in 1/U. We summarize 



the main results here. Details will be given elsewhere. 



In the SC channel, for small U r . one finds IT? 
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for arbitrary filling, with 
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7r Jo y 2p(/i) w c 



(8) 



and io c is the cutoff frequency for the square well (the second line holds when T c <C uo c ). 
In the limit uj c — > 0, one finds / — > \n(W/2u c ), as in the original work [|IJ, but the above 



expression also holds for arbitrary electronic DOS |18[ . 

Calculations in the CDW phase are more difficult. Restricting to the case of half filling 
(/i = 0) and again for small U c , one finds 



T. 



CDW 



T^ DW {0) 



exp 



exp 



1 



1 



p(0)[{2 - a(U c )}\U\ - U e ] p(0){2 - a(0)}|*7|J 
U c /\U\ + a(U c )-a(0) 



(9) 



{2-a(0)}M0)|tT| 

with a(U c ) a parameter that measures the reduction of the direct electron-phonon attraction 
by the exchange diagrams in Fig. 1(a). This parameter satisfies < a < 1 with a — > as 
Q — > and a — > 1 as Q — > oo. An estimate for a yields 

Ip(0)(\U\-U c 



a(U c 



1 



l-/p(0)|[/| 

which does approach 1 as c<j c — > oo (/ — > 0) and as u c 
Substituting Eq. (|10"D into Eq. (§) finally yields 

T C CDW (U C 



(10) 



[I ^l/p(0)(2\U\-U c )\ @ 
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exp 
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(ii) 



\U\^p(0)\U\ 

The resolution of the puzzle of how U c affects T C {CDW) versus T C (SC) is seen by ex- 
amining the small U c limits of Eqs. ((7|) and @. The pseudopotential effect disappears in 
the SC channel as U c — > 0, and the effect on T c is enhanced away from hal f filling since 



6 



p(/i) < p(0). On the other hand, the effect of U c on the CDW transition temperature is 
further reduced by the parameter / in Eq. flTT|). Thus, at least in the weak-coupling limit, 
U c initially reduces the SC T c more than it does the CDW T c . 

Fig. 3 shows the weak-coupling results (with the parameter / fitted to the QMC data for 
small Coulomb repulsion). One can see the weak-coupling formalism is excellent for g = 0.4, 
but becomes less accurate as the coupling strength increases. 

In the strong-coupling limit, where the electrons are paired into bipolarons at a temper- 
ature much higher than the transition temperature for the ordered phase, the initial effect 
of U c is to enhance the CDW T c , because Coulomb repulsion reduces the bipolaron binding 
energy Ef,, thereby increasing T c oc 1/Ef,. The analysis for the pure electron-phonon case, 
is easily modified by changing the energy of all intermediate states to take into account the 
Coulomb repulsion fi~9f. The CDW transition temperature at half filling then satisfies 



T^ DW {U C ) 1 

T CDW {0) « 1 _ ^ 



+ £(1 + S')(2 + S')...(n + S') 



1 



-sy 



S(l + 5)(2 + 5)...(n + 5)J (12) 
to second order in \U\, with S := |i7|/£2 and S' := S — (U c /Q). This analysis has been 

extended to fourth order, but it is too cumbersome to be shown here. The solid line in Fig. 3 

depicts the fourth-order approximation for T c (U c )/T c (0) [which lies with 10% of Eq. (|l|) for 

U c < 0.5\U\] when g — 1. Note how poor this approximation is as U c becomes a sizeable 

fraction of \U\. 

Conclusions. We have found that the conventional folklore for how Coulomb repulsion 
affects the electron-phonon interaction is flawed in assuming the SC transition temperature 
is more robust than the CDW transition temperature, but is resolved by detailed analysis 
of the weak-coupling theory. The Coulomb repulsion stabilizes the SC phase relative to 
the CDW phase and also allows for the appearance of incommensurate CDW order because 
T™ ax (SC) < T™ ax (ICDW) for finite U c . What happens to the phase diagram for U C >\U\1 
We conjecture that the CDW phase is taken over by the SDW phase at half filling, but do 
not know whether or not SC phases can remain stable away from half filling (due to the 



pseudopotential effect) or if paramagnetism prevails. Work along these lines is in progress. 
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FIG. 1. Lowest-order contribution to the irreducible vertex function in the (a) 
charge-density-wave, (b) superconducting, and (c) spin-density-wave channels. The solid lines 
denote electron propagators, the wiggly lines denote phonon propagators, and the dotted lines 
denote the Coulomb repulsion. 
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FIG. 2. Phase diagram for the Holstein-Hubbard model with g = 0.5, = 0.5, U = — 1.0, and 
U c = 0.0, 0.25, 0.5, and 0.75. The solid dots (lines) are for the commensurate CDW, the open dots 
(dotted lines) are for the incommensurate CDW, and the open triangles (dot-dashed lines) are for 
the SC as determined from a QMC (IPT) calculation. There are no transitions detected by the 
IPT for U c = 0.75. The dashed line through the QMC points is a guide to the eye. 
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FIG. 3. Relative change in the CDW transition temperature at half filling due to Coulomb 
repulsion. Four different cases are shown: g = 0.4 (solid dots); g = 0.5 (open dots); g = 0.625 (open 
triangles); and g = 1.0 (open squares). The dotted lines are the weak-coupling approximations for 
g = 0.4, 0.5, and 0.625. The parameter / is adjusted to produce the correct slope as U c — ► 0. The 
solid line is a fourth-order strong-coupling approximation for the g = 1.0 case. 
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